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Abstract. 

The Lagrangian properties of the velocity field in a magnetized fluid are studied 
using three-dimensional simulations of a helical magnetohydrodynamic dynamo. We 
compute the attracting and repelling Lagrangian coherent structures, which arc 
dynamic lines and surfaces in the velocity field that delineate particle transport in 
flows with chaotic streamlines and act as transport barriers. Two dynamo regimes are 
explored, one with a robust coherent mean magnetic field and one with intermittent 
bursts of magnetic energy. The Lagrangian coherent structures and the statistics of 
finite-time Lyapunov exponents indicate that the stirring/mixing properties of the 
velocity field decay as a linear function of the magnetic energy. The relevance of this 
study for the solar dynamo problem is discussed. 
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1. Introduction 

Transport in chaotic flows is governed by a combination of stirring and diffusion. 
Stirring refers to the transport, stretching, twisting and folding of fluid elements and, 
consequently, of scalar or vector quantities advected by the flow, such as temperature, 
light particles or magnetic field lines in a magnetohydrodynamic (MHD) system. This 
process creates complex tracer patterns in the flow, including filaments and sheets, 
as the fluid elements are deformed in different directions. Diffusion is responsible for 
homogenizing the distribution of tracers and blurring the patterns created by the chaotic 
stirring, being usually more important in small scales PQ. This paper deals with the 
problem of chaotic stirring in magnetized flows. Throughout the paper we use the terms 
"stirring" and "mixing" interchangeably. 

Passive scalars are quantities that are passively advected by the flow, i.e., their 
back-reaction on the advecting velocity field is disregarded. They constitute a powerful 
way to study transport in hydrodynamical and MHD flows (for a review, see Falkovich 
et al. [2]). We employ passive scalars to investigate how the magnetic field can 
affect particle transport and the stirring/mixing properties of a velocity field in MHD 
simulations through the Lorentz force. We adopt direct numerical simulations of resistive 
three-dimensional (3-D) compressible MHD equations with a helical forcing, which has 
been used elsewhere as a prototype of the a 2 dynamo model of mean field dynamo 
theory 014]. 

In the Lagrangian approach to turbulent transport the dynamics of fluids is studied 
by following the trajectories of a large number of fluid elements or tracer particles. 
The specific trajectories of individual particles are not very useful in this type of 
investigation in chaotic flows, since sensitivity to initial conditions means that particles 
that are arbitrarily close may experience exponential divergence with time. However, it 
is possible to detect certain material lines in the flow that repel or attract fluid elements. 
These repelling and attracting material lines are time-dependent analogous to the stable 
and unstable manifolds of hyperbolic fixed points in dynamical systems theory and form 
transport barriers in flows with chaotic streamlines, being called Lagrangian coherent 
structures (LCS). The LCS have been used to describe hydrodynamic turbulence in 3-D 
numerical simulations [5] , laboratory experiments [HI [7] and observational data of oceans 
[8j |9] and the atmosphere [10] , as well as 2-D numerical simulations of magnetized fusion 
plasmas [UJ, magnetic reconnection [12] , and 3-D MHD simulations of conservative JT3] 
and dissipative [4] fields. 

One of the most widely used Lagrangian tools are the finite-time Lyapunov 
exponents (FTLE), also known as direct Lyapunov exponents. The FTLE are a measure 
of local chaos and quantify the dispersion of particles in a region of the flow during a 
finite time. In the context of dynamo theory, the stretching rate of material lines in 
a fluid can be used to explain the amplification of magnetic fields by the mechanism 
of stretch-twist-fold dynamo [14). Examples of applications of the FTLE in dynamo 
simulations include the growth of seed magnetic fields in the kinematic dynamo problem 
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[T5l [T6l [T7j . nonlinear MHD dynamos [THl [HI [20] and the amplification of interstellar 
magnetic fields and turbulent mixing by supernova-driven turbulence in compressible 
MHD simulations [21 J . It was shown by Haller [22] that FTLE can also be used to 
identify repelling and attracting Lagrangian coherent structures. 

We present the detection of Lagrangian coherent structures for two different dynamo 
regimes in the 3-D compressible MHD equations with the isotropic and helical ABC 
forcing. We focus on the change in transport and mixing properties of the flow when 
the system undergoes a transition whereby a large-scale spatially coherent magnetic 
field loses its stability. The transition, which occurs after an increase in the magnetic 
diffusivity, results in strongly intermittent time series of magnetic energy. In the 
intermittent regime, the lower magnetic energy causes chaotic mixing to increase, 
resulting in higher stretching rates of material lines. Chaotic mixing is quantified by 
the finite-time Lyapunov exponents, which show a linear dependence on the magnetic 
energy. In section II of this paper we define Lagrangian coherent structures and how 
they relate to chaotic stirring in fluids; Section III describes the model of MHD dynamo 
adopted; the numerical analysis is discussed in section IV and section V presents the 
conclusions and possible ways to apply our techniques to observational data of the solar 
dynamo. 

2. Lagrangian Coherent Structures 

Let D C I 3 be the domain of the fluid to be studied, let x(t ) € D denote the position 
of a passive particle at time to and let u(x, t) be the velocity field defined on D. The 
motion of the particle is given by the solution of the initial value problem 

x = u(x,t), x(£ ) = x . (1) 

Let us define the following flow map: 4 > t° +T '■ x (^o) ^ x (^o + T )- The deformation 
gradient is given by J = d0^ +r (x)/dx and the finite-time right Cauchy-Green 
deformation tensor is given by A = J T J. Let Ai > A 2 > A 3 be the eigenvalues of 
A. Then, the finite-time Lyapunov exponents or direct Lyapunov exponents of the 
trajectory of the particle are defined as [23] : 

(7* 0+T (x) = In y/Xi, i = 1,2, 3. (2) 
|t| 

The maximum FTLE gives the finite-time average of the maximum rate of 
divergence or stretching between the trajectories of a fiducial particle at x and its 
neighboring particles. The maximum stretching is found when the neighboring particle 
y is such that 5x = x — y is initially aligned with the eigenvector of A associated 
with Ai. A positive a± is the signature of chaotic streamlines in the velocity field. The 
other exponents provide information about stretching/contraction in other directions 
and can be useful to interpret the local dynamics of the fluid. In an ideal conductive 
fluid, the frozen-in condition implies that a magnetic line aligned with an infinitesimal 
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vector connecting two close fluid elements will evolve as this vector [21] . As pointed 
out by Balsara et al. [21], for finite resistivity and compressible flows, flow regions with 
three positive Lyapunov exponents expand in all three directions and tend to dilute 
out the magnetic field; regions with two positive Lyapunov exponents and one negative 
exponent tend to concentrate the magnetic fields into sheet-like structures; regions 
with one positive and two negative exponents tend to mold the magnetic fields into 
filamentary structures; compression in all directions is found when all exponents are 
negative. On the other hand, local minima in the maximum FTLE field might provide 
a way to detect the position of the center of vortices in the velocity field, since vortices 
may be viewed as material tubes of low particle dispersion [25] . 

Finite-time Lyapunov exponents are also useful to detect attracting and repelling 
material lines that act as barriers to particle transport in the velocity field. A material 
line is a smooth curve of fluid particles advected by the velocity field [22]- These 
attracting and repelling material lines are the analogous of stable and unstable manifolds 
of time-independent fields. The study of 2-D flows is helpful to understand the role of 
material lines. Consider a 2-D steady flow, where the velocity field does not change 
with time. In the presence of counter-rotating vortices, hyperbolic (saddle) points are 
expected to be found, such as the one illustrated in Fig. [jja). The trajectories of 
passive scalars follow the velocity vectors in the vicinity of the hyperbolic point. Thus, 
particles lying on the stable manifold (green line) are attracted to the saddle point in the 
forward-time dynamics and trajectories on the unstable manifold (red line) converge to 
the saddle point in the backward-time dynamics. 

Two particles are said to straddle a manifold if the line segment connecting 
them crosses the manifold. The maximum FTLE has particularly high values on the 
stable manifold in forward-time, since nearby trajectories straddling the manifold will 
experience exponential divergence when they approach the saddle point, as shown in 
Fig. [ijb). Similarly, the FTLE field exhibits a local maximizing curve (ridge) along the 
unstable manifold in backward-time dynamics, since trajectories straddling the unstable 
manifold diverge exponentially when they approach the saddle point in reversed-time, 
as in Fig. [l](c). Thus, ridges in the forward-time FTLE field mark the stable manifolds 
of hyperbolic points and ridges in the backward-time FTLE field mark the unstable 
manifolds. 

Analogously, for a time-dependent velocity field, regions of maximum material 
stretching generate ridges in the FTLE field. Thus, repelling material lines (finite- 
time stable manifolds) produce ridges in the maximum FTLE field in the forward-time 
system and attracting material lines (finite-time unstable manifolds) produce ridges 
in the backward-time system [TT| |22| 123] . These material lines are called Lagrangian 
coherent structures (LCS). 

Stable and unstable manifolds are invariant sets, which means that a particle on 
the manifold will stay on it for all time. They form natural barriers to transport 
between different regions of a fluid, as seen in Figure [2] Figure |2](a) shows a topological 
configuration where one branch of the unstable manifold of the saddle point S\ smoothly 




Figure 1. (Color online) (a) Velocity vector field with a saddle point (black circle) 
and its stable (green) and unstable (red) manifolds; (b) schematic drawing of the 
forward-time trajectories (dashed lines) of two passive particles that straddle the stable 
manifold of a saddle point; (c) schematic drawing of the backward-time trajectories 
(dashed lines) of two passive particles that straddle the unstable manifold of a saddle 
point. 



joins the stable manifold of another saddle point S 2 in a heteroclinic connection. 
Simultaneously, the two branches of the unstable manifold of S 2 are connected to 
the stable manifold of Si, enclosing regions A and B. Particles trapped in A or B 
cannot cross the barriers formed by the manifolds, since these are invariant sets. Figure 
[2]^b) shows another type of trapping region, formed by a homoclinic connection, where 
one branch of the unstable manifold of a saddle point joins its own stable manifold. 
Trajectories in regions A and B usually circulate around a focus, as the manifolds mark 
the borders of vortices in the velocity field. Transport between different vortices is only 
possible when there is a transversal crossing between stable and unstable manifolds, 
through a mechanism called lobe dynamics [26, 27J. It is easier to understand this 
mechanism with a periodic flow. Suppose that the velocity field is time-dependent 
but periodic, such that u(x, t) = u(x, t + T), where T is the period. Let F be the 
stroboscopic Poincare map defined by F(x(t)) = 0j +T (x). There are still points where 
the velocity is instantly zero, but now they are moving. Since these points are not fixed, 
they are called stagnation points. After T time units a stagnation point will return to 
its original position. Therefore, under the map F a stagnation point is seen as a fixed 
point. Figure [2]^c) shows a heteroclinic tangle, where two hyperbolic fixed points of F 
have associated stable and unstable manifolds which intersect in a number of points, 
forming a set of lobes that protrude from one region to the other. At time t , lobes 
Ai and A 2 belong to region A and lobes Bi and B 2 to region B. Particles trapped 
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in each lobe cannot cross their bordering manifolds, but as time goes by the dynamic 
manifolds are transported and deformed by the flow, since they are material lines. After 
one period, lobe A\ is mapped onto the lobe marked as F(Ai), which belongs to region 
B. The same happens to lobe A2, which is mapped onto F{A2). Further iterations 
of the Poincare map F may cause lobes F(Ai) and F(A 2 ) to be deeply immersed into 
region B. Similarly, lobes B\ and B2 in region B are mapped to lobes F(Bi) and F(B2) 
in region A. 



a) b) 




Figure 2. (Color online) (a) Hctcroclinic connection between saddle points Si and 
S2. There can be no transport of scalars between regions A and B. (b) Homoclinic 
connections. Again, there is no transport between regions A and B. (c) Heteroclinic 
tangle between the unstable manifold of Si and the stable manifold of S2. The 
transport of particles between regions A and B is possible through lobe dynamics. 



3. The model 

A compressible isothermal gas is considered, with constant sound speed c s , constant 
dynamical viscosity /i, constant magnetic diffusivity rj, and constant magnetic 
permeability //q- The following set of compressible MHD equations is solved 



d t In p + u ■ V In p + V • u = 0, 

d t u + u-Vu = -Vp/p + J x B/p + (fi/p) (V 2 u + VV ■ u/3) + f , 
d t A = u x B — 7//i J, 



(3) 
(4) 
(5) 
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where p is the density, u is the fluid velocity, A is the magnetic vector potential, 
J = V x B//i is the current density, p is the pressure, f is an external forcing, and 
Vp/p = c^Vlnp, where c 2 s = •yp/p is assumed to be constant. Nondimensional units 
are adopted, such that c s = po = po = 1? where p = (p) is the spatial average of p. 
Equations @-@ are solved with the PENCIL CODE [J in a box with sides L = 2tt and 
periodic boundary conditions, so the smallest wavenumber is ki = 1. The time unit 
is {c s ki)~ l and the unit of viscosity v and magnetic diffusivity rj is c s /ki. The initial 
conditions are In p = u = 0, and A is a set of normally distributed, uncorrelated random 
numbers with zero mean and standard deviation equal to 10~ 3 . The forcing function f 
is given by the strongly helical ABC flow, 

f(x) = —4 [sin kfZ + cos kfy, sin kjx + cos kfZ, sin kfy + cos kfx], (6) 
V3 

where Af is the amplitude and kf the wavenumber of the forcing function. 

Following Rempel et al. [3111], we use Af = 0.1, kf = 5, and the numerical resolution 
varies between 64 3 and 128 3 . Spatial averages are denoted by (■) and time averages by 
(•) r References to kinetic (-Re) and magnetic (Rm) Reynolds numbers are based on the 
forcing scale 

Re = \fU/u, Rm = XfU/r], (7) 

where v = p/po is the average kinematic viscosity, Xf = 2ir/kf is the forcing spatial scale, 
and U = (u 2 ) 1 ^ 2 is the mean velocity at a time when the magnetic field is saturated. 
The turnover time r = \f/u rms varies between r w 3 and r w 4.5 for our range of rj. 

4. Results 

4-1. Bifurcation diagrams 

We choose rj as the control parameter and fix v = 0.005, which in the absence of magnetic 
fields corresponds to a spatiotemporally chaotic flow with Re ~ 100. Figure [3|a) shows 
the bifurcation diagrams for the time-averaged magnetic ((E m ) t , red triangles) and 
kinetic ((Ek) t , black circles) energies as a function of rj (lower axes) or Rm (upper 
axes). Averages are computed after an initial transient is dropped. For large values 
of r), the seed magnetic field decays rapidly and there is no dynamo. At the onset of 
dynamo action at rj ~ 0.053 (Rm ~ 9.5), the magnetic energy starts to grow at the 
expense of kinetic energy, until it saturates. Figure |3^b) shows in the upper panel the 
time-averaged kinetic helicity, if & = (u • u) , where u = V x u is the vorticity, and in 
the lower panel the time-averaged magnetic helicity, H m = (A ■ B). For helically forced 
flows, the magnetic helicity is expected to have the same sign as the kinetic helicity 
in scales smaller than the energy injection scale and the opposite sign in larger scales 
[28} . Most of the magnetic helicity in our simulations is concentrated in large scales, 
as happens with the magnetic energy [3], thus, H m has the opposite sign as Hk in Fig. 

\ http://pencil-code.googlecode.com 
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[3[b). These quantities are crucial for the emergence of a large-scale mean field, as they 
are related to the a-effect in mean-field dynamo theory, which is responsible for the 
generation of a mean electromotive force along the mean magnetic field by turbulent 
fluctuations of the velocity and magnetic fields [29], [30]. The presence of kinetic helicity 
is thought to be responsible for the inverse transfer of magnetic energy from small 
scales to large scales, as well as the inverse transfer of magnetic helicity from the energy 
injection scale to larger scales [28J. In Fig. ^b), H k is high for large values of r\ and H m 
is null, since there is no dynamo and a maximally helical (Beltrami) forcing is applied 
to the flow. After the onset of dynamo, the magnetic field starts to contribute to the 
flow dynamics through the Lorentz force (second term at the right in Eq. [4| and 
decreases with rj, as \H m \ grows. 

We focus on two values of rj. For rj = 0.01 the magnetic field is close to equipartition 
after saturation, as seen in the comparison between the time series of B Tms and w rms in 
Fig. [4] For t] = 0.05, close to the onset of dynamo action, the magnetic energy is almost 
an order of magnitude smaller than the kinetic energy and the time series are strongly 
intermittent, as shown in Fig. [5] Two pairs of vertical lines in Fig. [5]mark the beginning 
and apex of two bursts of magnetic energy around times t = 2400 and t = 6000. This 
is a type of on-off intermittency due to a blow-out bifurcation, as discussed by Rempel 
et al. [3]. 

The magnetic field structures are depicted in Fig. [6] for the two values of rj and 
different times. For rj = 0.01 (upper panel) there is a robust coherent large-scale B z 
component accompanied by small-scale turbulent fluctuations. For 77 = 0.05 (lower 
panel), the magnetic field displays an intermittent switching between coherent and 
incoherent large-scale structures [Fig. Knb)] and there is no preferred direction for field 
alignment. 



4-2. Lagrangian analysis 

The contrast between the Eulerian and Lagrangian analyses is depicted in Figure [7] for 
77 = 0.01, when the magnetic field has settled to a spatially regular mean field. Figure 
[7^ a) shows the line integral convolution (LIC) plot [3T] for the velocity field at z = 0. 
The LIC plot reveals the streamlines of the (u x , u y ) velocity components on this plane 
at t — 1000. Arrays of counter-rotating vortices found in the ABC-flow can still be 
seen, intermixed with long streaks and displaced vortices. To obtain the Lagrangian 
coherent structures we compute the maximum FTLE field. For the attracting LCS 
(finite-time unstable manifolds) we need to integrate the compressible MHD equations 
backward in time, which is a major problem, since the system is dissipative. We resort 
to interpolation of recorded data to compute these fields. A run of Eqs. Q-Q from 
to — t to to + T is conducted and full three-dimensional snapshots of the velocity fields 
are saved at each dt = 0.5 time interval. Following [11] and [32], linear interpolation in 
time and third order Hermite interpolation in space is used to obtain the continuous set 
of vector fields necessary to obtain the particle trajectories by integration of Eq. ([I]). 




Figure 3. (Color online) Bifurcation diagrams as a function of r\: (a) kinetic (black 
circles) and magnetic (red triangles) energies; (b) kinetic (black circles) and magnetic 
(red triangles) helicities. 



For backward integration, x = — u(x, t) is solved instead, as snapshots are read from to 
to to — r - Figure [7^b) shows the backward-time maximum FTLE field computed with 
t = — 10 and to — 1000. Bright colors correspond to large values of o\ and dark regions 
to low values. The ridges seen as bright red lines approximate the attracting LCS. Figure 
|7^c) shows the forward-time maximum FTLE field for r = 10, whose ridges provide the 
repelling LCS. Figure |7[d) is a superposition of (b) and (c), and represents the so- 
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Figure 4. (Color online) Time series for u rms (upper panel) and B rms (lower panel) 
for v = 0.005 and rj = 0.01. 
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Figure 5. (Color online) Time series for u rms (upper panel) and B lms (lower panel) 
for v = 0.005 and rj = 0.05. 



called "Lagrangian skeleton of turbulence" [32]. Figures [7^e) and[7^f) are enlargements 
of the square regions in Figs. [7^a) and[7](d), respectively. Note that the LIC plot of the 
velocity field in Fig. [7](e) shows a structure similar to the homoclinic connections of Fig. 
2]^b). Here, the arrows point to two hyperbolic stagnation points in the (u x ,u y ) field. 
On the other hand, when one moves to the Lagrangian frame (Fig. [7|f)) the picture 
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Figure 6. (Color online) (a) Intensity plots of B z at four different times, showing 
the evolution of a large-scale coherent pattern modulated along the x direction for 
r] = 0.01; (b) same as (a) but for r\ — 0.05, showing intermittent switching between 
ordered and disordered patterns. 

becomes much more complex, with a number of homoclinic and heteroclinic crossings, 
as in Fig. [2]^c). The two larger arrows point to the same location of the stagnation 
points. The smaller arrows indicate two lobes that cross other LCS and permit the 
transport of particles between vortices through lobe dynamics. The LCS in Fig. [7] were 
computed using 384 x 384 fiducial particles uniformly distributed on the plane z = 0. 
For each fiducial particle, the trajectories of six near neighboring particles are computed 
to obtain the deformation gradient by second-order centered finite-differences. 

In order to quantify the degree of particle dispersion or chaotic mixing in the flow, 
Fig. [8] shows the probability distribution functions (PDFs) of the three finite-time 
Lyapunov exponents for the same state shown in Fig. [7} The PDFs were obtained from 
a set o 64 3 initial conditions uniformly distributed in the box at to — 1000, with r = 10. 
There is a considerable amount of trajectories with two positive Lyapunov exponents, 
thus sheet-like magnetic field structures are expected. The broad tails in o\ are due to 
initial conditions that are very close to the repelling material lines, which are regions 
where the stretching is stronger than the average. On the other hand, the broad tails in 
negative values of 03 reflect contraction in the vicinity of the attracting material lines. 
Since the flow is weakly compressible, with the Mach number below 0.4, for almost all 
initial conditions one of the exponents is close to zero. The PDF for 02 shows a Gaussian 
distribution. The overbars on denote average values. 

For 77 = 0.05 the time series of B ims and u rms are intermittent. To understand the 
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e) 








Figure 7. (Color online) (a) Line integral convolution plot showing the streamlines 
of the ccy-components of the velocity field at t = 1000 for rj — 0.01; (b) the 
attracting Lagrangian coherent structures (LCS) (red); (c) the repelling LCS (green); 
(d) superposition of (b) and (c); (e) enlargement of the square region of (a); (f) 
enlargement of the square region of (d). 
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Figure 8. PDFs of the finite-time Lyapunov exponents at to — 1000 and r = 10 for 
T] = 0.01. 



influence of B on u, the FTLE are computed for different initial times marked with 
vertical lines in Fig. |5j Figure [9] shows the LIC and LCS plots at t — 2100, just before 
a burst of magnetic energy in the time series of Fig. [5Fb). In comparison to rj — 0.01, 
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there seems to be less order in the distribution of vortices in the LIC plot of Fig. |9[a) 
than in [71(a) and the greater complexity in the distribution of material lines in the LCS 
plots of Figs. |9](b)-(d) and (f) indicates that transport of passive scalars is enhanced 
due to the frequent crossings of attracting and repelling lines. This is expected, since 
-B rms is much lower for 77 = 0.05 and has a smaller impact on the velocity field. 



r) = 0.05 t = 2100 




Figure 9. (a) Line integral convolution plot showing the streamlines of the xy- 
components of the velocity field at t = 2100 for 77 = 0.05; (b) the attracting Lagrangian 
coherent structures (LCS) (red); (c) the repelling LCS (green); (d) superposition of 
(b) and (c); (e) enlargement of the square region of (a); (f) enlargement of the square 
region of (d). 



At t = 2400 the time series of B Tms has a peak of energy burst. As seen in Fig. 
10, there is a stronger impact of this magnetic energy release on the velocity field in 

The LIC plot of Fig. 10 a) does not show much 
However, the LCS plots of Figs. 10 ^b)-(d) show 



comparison with t=2100 (Fig. [9 
difference in relation to Fig. |9j(a 



wider regions of low particle dispersion. This is clearer in Fig. 10 f), which shows an 
intermediate level of complexity in comparison to Figs. [7](f) and[9]T). Thus, a stronger 
magnetic field diminishes chaotic mixing in the velocity field, which is measured by the 



PDFs of the FTLE, shown in Fig. 11 



One can see that the PDFs for the intermittent dynamo (Fig. |lT[a),(b)) are wider 
than for the regular mean-field dynamo (Fig. [8]). They also have a larger a±, revealing 
that the flow is more chaotic for 7/ = 0.05 than for rj = 0.01. Moreover, broader tails in 
the PDFs of Fig. 11 mean that more intermittency is to be expected in the evolution 
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of passive scalars at rj = 0.05. A summary of the results is found in Table [TJ which 
shows o"i,2,3 and their standard deviations for rj = 0.01 at t — 1000 and for r] = 0.05 at 
four values of t representing the beginning and apex of the two magnetic energy bursts 
indicated in Fig. [5j The mean value o\ and the standard deviation decrease at both 
bursts. Figure [12] is a plot of <7i as a function of -B rms using only data from Table [TJ 
which are fitted with the linear equation o\ ~ 0.348 — 0.345-B rms . Although we need 
more statistics to draw conclusions, our preliminary results suggest that the decay of 
chaoticity in the velocity field is proportional to B ims . 



r| = 0.05 t = 2400 




X 



Figure 10. (a) Line integral convolution plot showing the streamlines of the xy- 
components of the velocity field at t — 2400 for rj = 0.05; (b) the attracting Lagrangian 
coherent structures (LCS) (red); (c) the repelling LCS (green); (d) superposition of 
(b) and (c); (e) enlargement of the square region of (a); (f) enlargement of the square 
region of (d). 



5. Conclusions 

We have employed Lagrangian coherent structures (LCS) and the statistics of finite-time 
Lyapunov exponents (FTLE) to study chaotic stirring in 3-D MHD dynamo simulations 
with helical forcing. Attracting LCS provide the pathways that are more likely to be 
followed by passive scalars and their crossings with repelling LCS provide the mechanism 
for transport between different regions of the fluid. The PDFs of FTLE provide a 
quantification of chaotic mixing in the flow. We explored the impact of the magnetic field 
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Figure 11. PDFs of the finite-time Lyapunov exponents for 77 = 0.05 at (a) i = 2100 
and (b) t = 2400. 



Table 1. Mean values and standard deviations of finite-time Lyapunov exponents. 





r] = 0.01 


r] = 0.05 




t = 1000 


t = 2100 


t = 2400 


* = 5600 


* = 6000 




0.249 


0.328 


0.298 


0.332 


0.303 


0-2 


0.006 


0.030 


0.017 


0.032 


0.021 


0-3 


-0.245 


-0.328 


-0.292 


-0.333 


-0.299 


std(a"i) 


0.096 


0.125 


0.123 


0.125 


0.124 


std(a 2 ) 


0.067 


0.098 


0.089 


0.098 


0.091 


std(cr 3 ) 


0.098 


0.137 


0.124 


0.138 


0.127 



on the velocity field in a saturated nonlinear dynamo and in an intermittent dynamo, 
and the maximum FTLE was shown to be a linear function of the magnetic energy. 
The increase in the flow's chaoticity when the magnetic diffusivity is increased from 
i] = 0.01 to i] = 0.05 is the result of the reduction of the effect of the Lorentz force 
upon the velocity field. Enhanced chaoticity leads to stronger line stretching and field 
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Figure 12. Maximum finite-time Lyapunov exponent of the velocity field (01) as a 
function of B rms using data from Table [T] The fitted line has slope 7 ~ —0.345. 

amplification, and the "competition" between this effect and destruction of magnetic 
flux due to magnetic diffusion seems to be the main cause of the intermittent time series 
of magnetic energy observed when rj = 0.05, which is close to the critical value for 
dynamo action. 

Our analysis has direct applications in astrophysics, where the equipartition- 
strength magnetic fields observed in planets and stars are thought to be the result 
of a dynamo process, whereby kinetic energy from the motion of a conducting fluid is 
converted into magnetic energy [33]. Experimental detection of LCS and computation 
of FTLE in the solar surface can be performed using velocity fields estimated from 
observational data. Such estimations can be obtained from digital images using the 
optical flow algorithm, employed by [34J to extract the velocity field from images of 
coronal mass ejections obtained with the SOHO LASCO C2 coronagraph. Recently, 
horizontal velocity fields in the photosphere were inferred from Hinode images [33] 
and the Swedish Vacuum Solar Telescope [36]. Solar subsurface flows can be inferred 
from helioseismic data [37], thus LCS can also aid the tracing of particle transport by 
turbulence in stellar interiors. 
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